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1 Introduction 

Let p be a rational polytope in V = M. d and h(x) a polynomial function on 
V. A classical problem in Integer Programming is to compute the sum of 
values of h(x) over the set of integral points of p, 

s(p,h)= h ^)- 



When p is dilated by an integer n e N, we obtain a function of n which is 
quasi-polynomial, the so-called Ehrhart quasi-polynomial of the pair (p, h) 



d+N 



S(np, h) = ^ E m (p, h, n)n" 

m=0 



of degree d + N where iV = degh. The coefficients E m (p, h,n) are periodic 
functions of n G N, with period the smallest integer q such that qp is a lattice 
polytope. 

Replacing h(x) by an exponential, we are led to study the analytic func- 
tion on V* defined by 

s(p)(o= E eiM - 

xepnz d 

If p is any rational polyhedron, this sum still makes sense as a meromorphic 
function defined near and the map p i— ► S'(p)(^) is a valuation. 
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In [1], we proved that the meromorphic function 5(p)(^) has a local Euler- 
Maclaurin expansion 

s(pxo = I>(t(M))(o l e ^ x)dx - 

The sum is taken over the set of faces f of the polyhedron p. For each 
face f, the function //(t(p, f))(£) is holomorphic near 0, and it depends only on 
the transverse cone t(p , f ) of p along f . More precisely, once a rational scalar 
product is chosen on V, we define canonically a map a i— > /i(a) from the set 
of rational affine cones a in quotient spaces of V, with values in the space of 
functions on V which are analytic near 0, then we prove that these functions 
satisfy the above formula. The map a i— > /i(o) is invariant under lattice 
translations, equivariant with respect to lattice preserving isometries, and it 
is a valuation on the set of affine cones with a fixed vertex ([1], Theorems 17 
and 18). 

It is easy to see that the Ehrhart quasi-polynomial can be computed in 
terms of the Taylor coefficients of the functions /u(t(p, f))(£)- For example, if 
p is a lattice polytope and h(x) = 1, we have ([3], Corollary 28) 

(1) Card(npnZ d ) = ^ /U (t(p,f))(0)vol(f)n dim f. 

f 

Using the valuation property of /i(ct) and Barvinok's decomposition of a 
cone into unimodular cones, we thus obtained in [4j an algorithm for com- 
puting the Ehrhart quasi-polynomial. It has polynomial length with respect 
to the input (p, h), when the dimension d and the degree N are fixed. 

The valuation S(p, h) has a natural generalization used by Barvinok in 
[2], the mixed valuation S L (p, h), where L C V is a rational vector subspace. 
Denote the projected lattice on V/L by A^. For a polytope p C ^ and a 
polynomial h(x) 

S L (p,h) = I Kx)dx. 

yeAv/L Jpn(y+L) 

In other words, the polytope p is sliced along lattice affine subspaces parallel 
to L and the integrals of h over the slices are added up. For L — V, there 
is only one term and ^(p, h) is just the integral of h(x) over p, while, for 
L = {0}, we recover S(p, h), the sum of values of h(x) over the set of integral 
points of p. 
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In the case h(x) = 1, we write S(p) and S L (p) in place of S(p, 1) and 
S L (PA). 

Using these mixed valuations, Barvinok gave an algorithm which com- 
putes the r + 1 highest degree Ehrhart coefficients of S(np) = Card(npflZ d ), 
when p is a simplex in IR d . Barvinok' s algorithm has polynomial length when 
d is an input, provided r is fixed. The method consists in reducing the 
problem to summations over lattice points in dimension < r. 

Barvinok considers particular linear combinations 

Lec 

where £ is a finite set of rational vector subspaces of V which is closed under 
sum, and the coefficients p(L) are integers which satisfy the following relation 
between characteristic functions: 

Lec 

where L 1 - C V* is the orthogonal of L. We call a function L — > Z with this 
property a patchwork function on £. 

When p is dilated by an integer n, S L (np) is again given by a quasi- 
polynomial in n, as is a linear combination 

d 

Y,P(L)S L (np) = J2"m(P,n)n m . 

LeC m=0 

The main theoretical result of [2], Theorem (1.3), is the following: if L 
is a family of subspaces which is closed under sum and contains the vector 
subspace lin(f) parallel to f, for every face f of codimension < r of p, and if p is 
a patchwork function on C, then the r + 1 highest degree coefficients v m {p-, n), 
for m — d, . . . , d — r, are equal to the corresponding Ehrhart coefficients of 
S(np) = Card(np nZ d ). 

In the present article, we introduce the meromorphic functions which 
extend S L (p). For any polyhedron p, 

s L (pm= E / e{Mdx 

yeAv/L Jpn(y+L) 
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is defined as a meromorphic function near £ = 0. We show that S' L (p)(£) 
also enjoys a local Euler-Maclaurin expansion (Theorem |SJ), 

s L (p)(o = 5> L wp'f))(£) [ e ^ x)dx - 

Furthermore, for a linear combination of Barvinok type, if f is a face 
of p such that lin(f) e £, we prove that the f-term in the Euler-Maclaurin 
expansions of 

^G>)(0 = £p(^(P)(£) 
and of the tzsna/ valuation S(p)(£) are equal (Theorem H7|) : 
J]p(L)^(t(p,f))(0=/i(t(PJ))(0- 

This is the main result of the present article. From the relation between 
Ehrhart quasi-polynomials and Euler-Maclaurin expansions, it implies Barvi- 
nok's Theorem (1.3). 

Actually, we derive from Theorem [T7] another computation of the r + 1 
highest coefficients of the Ehrhart quasi-polynomial for a pair (p,h), based 
on Brion's decomposition of a polytope into cones, in the line of [1] and [6]. 

For each vertex s of p, let c s be the cone of feasible directions of p at 
s. Instead of the full family C generated by taking sums of the subspaces 
lin(f), when f runs over the set of faces of codimension < r of the polytope p, 
we consider, for each vertex s of p, the family C s generated by faces of c s of 
codimension < r. The point in taking a family which depends on s lies in the 
case where p is simplicial. Then C s consists only of the spaces lin(f) where 
f is a face of c s of codimension < r, as this set is already closed under sum. 
Moreover the coefficients p(L) are just signed binomial coefficients (Lemma 
[T5|) . and the computation of S L (c s ) is easier when L is parallel to a face of 
c s (Example [6]). 

Let us describe our method in the simpler case of a lattice polytope p and 
polynomial h(x) = 1. By Brion's theorem, we have 

s(p)(o = J2 e{Ms ^(0- 

s 

For each vertex s, let p s : C s — * Z be a patchwork function. We define 

B(p)(o = J2 e{MsCs,p8 ( c ^- 



4 



For the dilated polytope np, we have 

s m>0 s 

Hence, the meromorphic function ^XL(£> s ) m ^(CsXO * s actually regular 
at £ = and its value at £ = is the mth Ehrhart coefficient of p. 
We have similarly 

m>0 s 

The meromorphic functions 5* L (c s )(£) and 5* Ps (c s )(£) have a special form: 
they can be written as the quotient of an analytic function by a product 
of d' < d linear forms. Such a function <p has an expansion into rational 
functions (p = ^2j>-d^\j] wnere 4>\j\ is homogeneous of total degree j. 
Now it follows from our main theorem that, for m > d — r, we have 

S(c s ) hm] (0 = S c ^(c s ) hm] (Oi 

hence the zero degree terms of 5^ s (£, s) m S(t s )(g) and 5^ s (£, s) m S Cs,f)B (c s )(^) 
are equal. Therefore the latter is also analytic at £ = and its value at £ = 
is the mth Ehrhart coefficient of p. This is the content of Theorem [201 

Thus, besides taking care of any polynomial h(x), not only h(x) = 1, this 
method to compute the r + 1 highest degree Ehrhart coefficients for the pair 
(p, h) leads to a simpler algorithm than the one proposed in [2]. When p is a 
rational simplex, the contributions of the terms of the form S L (c s )(^) when 
L G C s are immediately reduced to the computation of a function S(a) with 
a a simplicial cone of dimension smaller or equal to r. 

There is also another possible implementation of an algorithm to compute 
the r + 1 Ehrhart highest degree coefficients for the pair (p, h) based on the 
results of [4J. As seen from Equation ((Tj), this involves the computation of the 
analytic function /i(t(p,f)), also associated to simplicial cones in dimension 
smaller or equal to r. We plan to compare the implementation of both 
methods in the near future. 
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2 Local Euler-Maclaurin expansion of a mixed 
valuation S L 



We consider a rational vector space V of dimension d, that is to say a finite 
dimensional real vector space with a lattice denoted by Ay or simply A. 
We will need to consider subspaces and quotient spaces of V, this is why 
we cannot just let V = R d and A = Z d . By lattice, we mean a discrete 
additive subgroup of V which generates V as a vector space. Hence, a lattice 
is generated by a basis of the vector space V. A basis of V which is a Z-basis 
of Ay is called an integral basis. The elements of A are called integral. An 
element x G V is called rational if qx G A for some integer q ^ 0. The space 
of rational points in V is denoted by Vq. A subspace L of V is called rational 
if L H A is a lattice in L. If L is a rational subspace, the image of A in V/L 
is a lattice in V/L, so that V/L is a rational vector space. We will call the 
image of A in V/L the projected lattice. 

Example 1 Let V = IR 2 standard lattice Z 2 . Lei i>i,i>2 6e £tt>o primitive 
integral independent vectors. Using an integral basis with first basis vector 
Vi, a straightforward computation shows that the projected lattice on IR 2 /!^ 
is Z det ^ v ^ , where v 2 is the projection of v 2 on R 2 /Ri>i. 

A rational space V, with lattice A, has a canonical Lebesgue measure, 
for which V/A has measure 1. An affine subspace L of V is called rational if 
it is a translate of a rational subspace by a rational element. It is similarly 
provided with a canonical Lebesgue measure. We will denote this measure 
by dm L . 

We will denote elements of V by latin letters x,y,v, . . . and elements of 
V* by greek letters £, a, We denote the duality bracket by (£, x). 

If S is a subset of V, we denote by < S > the affine subspace generated 
by S. If S consists of rational points, then < S > is rational. Remark 
that < S > may contain no integral point. We denote by lin(S') the vector 
subspace of V parallel to < S >. 

If S is a subset of V, we denote by S 1 - the subspace of V* orthogonal to 

S: 

S ± = {£ e V* ; (£, x) = for all x G 5}. 

If L is a subspace of V, the dual space (V/L)* is canonically identified with 
the subspace L L <zV* . 
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A convex rational polyhedron p in V (we will simply say polyhedron) is, 
by definition, the intersection of a finite number of half spaces bounded by 
rational afline hyperplanes. We say that p is solid (in V) if < p >= V. A 
polytope p is a compact polyhedron. 

The set of non negative real numbers is denoted by R+. A convex rational 
cone c in V is a closed convex cone Y^l=i ^+ v i which is generated by a finite 
number of elements i>j of Vq. In this article, we simply say cone instead of 
convex rational cone. 

An afline (rational) cone a is, by definition, the translate of a cone in V 
by a rational element s G Vq. This cone is uniquely defined by o. 

A cone c is called simplicial if it is generated by independent elements of 
Vq. A simplicial cone c is called unimodular if it is generated by indepen- 
dent integral vectors V\, . . . ,Vk such that {v\, . . . ,Vk} can be completed in 
an integral basis of V. An afline cone a is called simplicial (resp. simplicial 
unimodular) if it is the translate of a simplicial (resp. simplicial unimodular) 
cone. 

An afline cone a is called pointed if it does not contain any straight line. 

The set of faces of an afline cone a is denoted by jF(a). If a is pointed, 
then the vertex of a is the unique face of dimension 0, while o is the unique 
face of maximal dimension dim a. 

Let us recall the definition of the transverse cone t(p, f) of a polyhedron 
p along one of its faces f. Let x be a point in the relative interior of f. The 
cone of feasible directions of p at x is the set c(p,f) := {v G V ; x + ev G 
p for e > small enough}. It does not depend on the choice of x. We denote 
the projection V — > V/ lin(f) by 7r f . Then t(p, f) is the image 7r f (f + c(p, f)) 
of the afline cone f + c(p,f) in V/lin(f). It is a solid pointed afRne cone in 
the quotient space V/ lin(f) with vertex vr f (< f >). In particular, if v is a 
vertex of p, the transverse cone t(p,v) coincides with the supporting cone 
v + c(p,v) C V. 

If a is an afline cone and f is a face of a, then c(a, f ) = a + lin(f) and the 
transverse cone t(a, f) of a along f is just the projection 7r f (a) of a on V/ lin(f). 

Definition 2 Denote byTiiV*) the ring of analytic functions around G V*. 
Denote by AA(V*) the ring of meromorphic functions defined around G 
V* and by M. £ {V*) C M.(V*) the subring consisting of those meromorphic 
functions 0(£) such that there exists a product of linear forms D{^) with 

D(0(l>(0 G H(V*). 
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A function 0(£) G A4e{V*) has a unique expansion into homogeneous rational 
functions 

m>>— oo 

where m is the total degree. 

If P is a homogeneous polynomial on V* of degree p, and D a product 
of r linear forms, then is an element in A4i(V*) homogeneous of degree 
m = p — r. 

Let us recall the definition of the function /(p) G M-t(V*) associated to 
a polyhedron p, (see for instance the survey [3]). 

Proposition 3 There exists a map I which to every polyhedron p C V asso- 
ciates a meromorphic function with rational coefficients I(p) G Aii(V*), so 
that the following properties hold: 

(a) If p contains a straight line, then l(p)=0. 

(b) If £ G V* is such that |e^'^| is integrable overp, then 

I(P)(0 = [ e^dm <p> (x). 
Jp 

(c) For every point s G Vq, we have 

/( S + p)(0 = e^/(p)(0. 

(d) The map I is a simple valuation: if the characteristic functions x{Pi) 
of a family of polyhedra pi satisfy a linear relation ^j^x(pi) = ; then the 
functions I (pi) satisfy the relation 

{i,<Pi>=V} 

In the following proposition, we define the mixed valuation p i— > 5* L (p) 
associated to a rational vector subspace L Q V. To any polyhedron p, we 
associate a meromorphic function S L (p)(^) G M.{V*). If p is compact, this 
function is actually regular at 0, and its value for £ = is the valuation 
E L ±(p) considered by Barvinok [2]. 

We denote by A v / L the projection on V/L of the lattice A. 
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Proposition 4 Let L Q V be a rational subspace. There exists a map S L 
which to every rational polyhedron p C V associates a meromorphic function 
with rational coefficients S L (p) G Ai(V*) so that the following properties 
hold: 



(a) If p contains a line, then S L (p)=0. 




for every £ G V* such that the above sum converges. 

(c) For every point s G A, we have 

s L {s+pm = e^s L {pm. 

(d) The map S L is a valuation: if the characteristic functions x(pi) °f 
a family of polyhedra pi satisfy a linear relation ^\ Tix{pi) = 0, then the 
functions S L (pi) satisfy the same relation 

J2r l S L (p l ) = 0. 

i 

For L = {0}, we recover the valuation S given by 

S(P)(0 = £ e«-'\ 
xepnA 

provided this sum is convergent. 

For L = V, we have S v (p) — I(p), if p is solid, and S v (p) =0 otherwise. 

The proof is entirely analogous to the case L = {0}, see Theorem 3.1 in 
[3], and we omit it. 

Remark 5 The function S L (p) is actually an element of M. t {V*), but we 
do not prove it at this point. Let a be an affine cone and {vi} the generators 
of its edges. It will follow from the Euler-Maclaurin expansion of S L (a) 
(Theorem^) that Yli(^ v i)^ L ( a )(0 ^ s analytic near zero for any L. It would 
be interesting to prove it a priori. By Brion's theorem and the valuation 
property, it follows in particular that S L (p) G AiiiV*). 



9 



Example 6 Let a be a simplicial affine cone in the space V , and assume 
that L = lin(fi) for some face f 1 of a. In this case, S L (a)(C,) decomposes as 
product of an integral and a discrete sum. For simplicity, assume that a is 
solid. Let f 2 be the face of a such V = lin(fi) © lin(f 2 ). We write x = x\ + x 2 
and £ = £1 + £2 for the corresponding decompositions of x G V and £ G V* . 
Then a = cii + a 2 where cij is a simplicial affine cone in lin(fj). Let us 
denote by A 2 the projected lattice in V/ lin(fi) ~ lin(f 2 ) . From (dp, we obtain 
immediately 

(3) S L (a)fe + 6) = /(ai)(£i) e ^ X2) - 

Notice that the lattice A 2 is usually bigger than A D lin(f 2 ). 

Example 7 Let V = M 2 with the standard lattice. We compute S L (a) when 
a is a cone and L is a line. Let a = IR + t>i + IR + t>2, where v\, v 2 are two linearly 
independent primitive integral vectors. 

(a) Assume that L is the line supporting an edge of a, say L = M.vi. We 
identify V/L to M.v 2 . The projected lattice is A 2 = Z det( .^ v ^ , (Example^), 
hence, by §3±) in Example® we have 

(4) S L (a)(0- 



(e,"2> 

g I dct(v 1 ,v 2 )\ 



(b) Assume now that L is transverse to both edges of a. Assume that 
det(t>i, v 2 ) > 0. Let L = Ru where u is a primitive integral vector chosen so 
that det(M,f 2 ) > 0. Let cij = M. + u + WL + Vi for i — 1,2. We decompose the 
characteristic function of the cone a as x( a ) — x( a 2) + x( a i) — or 
x(<x) = ^(02) — x(ai) + i), depending on whether u belongs to a or not. 

Using the valuation property, case (a) and the relation 



1 



we obtain in both cases 
(5) S L (a)(0 = 



\S»J "7 V Y gdct(u,u 2 ) 



g dct(u,v-^ ) 



In this example, we see that {^,vi){^,v 2 )S L (a)(C,) is indeed analytic near 

e = o. 
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In the following theorem and its applications, we will consider the func- 
tions S L (p) when the space V is replaced with a quotient space W. We 
denote by C{W) the set of affine cones in W. Thus if a € C(W), and L a ra- 
tional subspace of W, the function S L (a) is a meromorphic function on W*. 
We are going to show that the function S L (a) has a local Euler-Maclaurin 
expansion, which generalizes the case L = {0} of [3]. 

Theorem 8 Let V be a rational space and Q a rational scalar product on 
V* . There exists a unique family of maps n^, indexed by pairs (W, L) where 
W is a rational quotient space of V and L is a rational vector subspace of W 
such that the family enjoys the following properties: 

(a) maps C(W) to H(W*), the space of analytic functions on W* . 

( b ) IfW = {0}, then /$>({()}) = 1. 

(c) For &mW > and L = W , then ^(a) = 0. 

(d) If the affine cone a G C(W) contains a straight line, then /i^(a) = 0. 

(e) For any affine cone a in W, the following formula holds 

s L (a)= &/£$ M H<*,mf) 

where the sum is over all faces of the cone a. 

In this last formula, the function /i^t 1 ' ^(t(a, f)) is considered as a 
function on W* itself by means of the orthogonal projection W* — > (W/ lin(f))* = 
(lin(f))- 1 - with respect to the scalar product on W* C V* . 

Proof. The proof is entirely similar to the case L = {0} studied in [1]. 
Note that [iff coincides with the map denoted by fj,w in [3J- The only new 
item is (c). It follows immediately from the relation S w (a) = 1(a). 

Remark 9 Let a be a solid cone in W , and let f be a face of a such that 
dimf < dimly. If L is transverse to the face f, that is, if L + lin(f) = W , 
then fJ^ff^ Mf) {t{a, f)) = 0. This follows from (c). 

From now on we omit the subscript W, thus we write \i L in place of /i^,. 
The next theorem and its proof are also entirely similar to the case L = {0} 
in g]. 
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Theorem 10 The analytic functions defined in Theorem^ have the follow- 
ing properties: 

(a) For any x G A, one has fi L (x + a) = /i L (a). 

(b) The map (a, L) \— > /i L (a) is equivariant with respect to lattice-preserving 
isometries. In other words, let g be an isometry of W which preserves the 
lattice A w . Then ^^(^(a))^ -1 ^) = . L ( &)(0- 

(c) If W is an orthogonal sum W = W\ © W 2 of rational spaces, Li C Wi 
and Oj is an affine cone in Wi for i = 1,2, then 

fi L ^(a 1 + a 2 )=^(a 1 )^(a 2 ). 

(d) For a fixed s G Wq, the map c — > fi L (s + c) is a valuation on the set 
of cones in W . 

(e) Let p C W be a rational polyhedron, then 

(6) S L (p)= /i L+lm{f)/lin(f) (t(PJ))/(f)- 

Example 11 Let us compute the function fi L for the various transverse 
cones of Example [7| We define a function B(u) on C, holomorphic near 
0, by 

1 1 

1 — e u u 

We have 

\&et(v 1 ,v 2 )\ 

I[m) = <£,«*>«,«*>■ 
Consider case (b) where L is transverse to both edges fi = IR+fj of a. 
Using the equation det(vi,v 2 )u = det(u, v 2 )vi — det(u,V\)v 2 , we have 

1 f det(u,v 2 ) det(«,vi)' 

/(0)({) = \ 

Thus we can rewrite (EP as 

S L (a) = //(a) + V L+rm(b)/rm{h \(t(a,U)m) +/(o), 



i=l,2 



with 

(7) A l (a)«) 



L \det(u, fi)/ \det(it, 
^x+iin (f i )/iiu(ft)((t( a> f •)) = o for i = 1, 2, 
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Observe that ([?]) is indeed regular at £ = 0. 

In case (a) where L = Ri> 1; we /iai>e L + lin(f 1 ) / lin(f 1 ) = {0}. Let us 
assume that det(t>i,t> 2 ) > 0. Then we have, by fcflj, 

V det(v 1 ,v 2 ) 

W ith a = #^4. 

As J(fi)(£) = ~ 7g^\ ; we can rewrite as 

S L (a) = ^(a) +/x{°>(t(a,f 1 ))/(f 1 ) + /i i+lin *)/ lin *)((t(a, f 2 ))/(f 2 ) + 1(a), 
with 



/i L (a)(0 ' 



det (v\ , v 2 ) / V det (wi , f 2 



which is indeed regular at £ = 0, and, again, 

M L+Bn(fe)/UB(fe)(( t ( aj f 2 )) =0. 

3 Barvinok valuations 

Let £ be a finite family of rational vector subspaces of V, and let p(L), Lg£, 
be a set of rational coefficients. The linear combination Els£ P(L)S L (p) is 
again a valuation on the set of polyhedra p C V, with values in _A/^(V*). By 
taking linear combinations, we obtain a local Euler-Maclaurin expansion for 
the function ElscPW^XO- 

Definition 12 Lei p C V be a polyhedron, 
(a) We denote 

Lec 

(b) We define the f-term in the local Euler-Maclaurin expansion of S^ c,p \p) 
to be 

^ p \t(p,f)) = ^p £ (^)^ +lm{f)/lin(f) (t(P,f)). 
Lec 
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Thus we have 



S <M(p)= /i (£ ' p) (t(PJ))/(f)- 



fe-^(p) 



We are going to compute the f-term in the case of the following particular 
linear combinations introduced by Barvinok [2]. 

Definition 13 The valuation S^ c ' p ^ is called a Barvinok valuation if 

(a) the family of sub spaces £ is stable under sum, 

(b) p is an integer valued function on the set £ such that the characteristic 
function of the union of the subspaces L 1 - Q V* can be written as the linear 
combination 



Definition 14 We call a function £ — > Z which satisfies (Q|) a patchwork 
function on £. 

As the set of orthogonal subspaces L L Q V* is stable under intersection, 
a particular function pc with this property can be computed in terms of the 
Moebius function of the partially ordered set £ ([7], vol I, section 3.7), as 
explained in [2j. 

Let us compute a patchwork function p in the following case. V = M. d with 
standard basis e^, i = 1, . . . , d, and £^ g is the set of subspaces Lj = ©jg/IRej 
with cardinal \I\ > q. The function p^ q defined below is actually the one 
associated to the Moebius function, but we will not need this fact. 

We denote the binomial coefficient , ,, m! M , by Ct,. 

k\(m—k)\ J rn 

Lemma 15 The function p^ q on £^ q defined by 



satisfies Equation 

Proof. If e % is the dual basis, the orthogonal space Lf is equal to 
^jRe'. Let £ = Ylt=i e ^LeCqL^- Let I be the set of indices 
i G [l,...,d] such that £j = 0. Then \I \ > q, and £ G Lj~ if and only 
I C J Q . Let \Iq\ = N. The value at £ of the right-hand side of (jE} is equal to 





Pd, q (Li) = {-ly-'ctx 



if \I\ = n, 



N 



E(N, q ) = y.p*m = 




ICIq n=q 
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We want to prove that E(N, q) = 1. Writing n = q + i, we have 

N-q 

E(N,q) = J2(-n i 



i=0 



(q + i)(N - q - i)W.(q - 1)!' 



Let us compute (g — 1)! (e(N + 1, q) — E(N, g)J . This is equal to 
1 J (JV + l)(jV+l_ g )| ^ + + 1 -g-i)! (iV-g-i)! 



i=0 

N+l-q 

r_ iv 

z!(JV + l-g-i)! 
We obtain E(JV, q) = E{q,q) = 1. □ 



EiV> 
r— iv — = o. 
1 ^ MAT -1-1 

i=0 



In the case of a Barvinok valuation, it turns out that the f-term in the 
Euler-Maclaurin expansion of S^ c ' p \p) coincides with that of S(p), if the 
vector subspace lin(f) belongs to the set C This is the crucial result of 
the present article. It is an easy consequence of the following combinatorial 
lemma. 

Lemma 16 Let C be a finite family of vector subspaces of V , stable under 
sum and let p be a patchwork function on C. 

(a) Let L 6 C. Then 

£ P(L) = 1. 

{L(LC, LCLq} 

(b) Let Lq ^ L\ be two subspaces in the family C. Then 

£ P{L) = 0. 

{LeC, L+L =L 1 ] 

Proof. There exists a £o £ Lq such that, for L E C, C,q E L ± if and only 
if L C Lq. We obtain (a) by evaluating both sides of (jSJ) at this particular 
element £o- Next, we deduce (b) from (a), by induction on dimLi — dimL . 
For two subspaces M C M' in the family C, let us denote 

f(M,M')= Yl P(L)- 

{LeC, L+M=M'} 
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If M = M>, we have f(M, M) = E {ieA zcM} P(L) = 1 by (a). 
We apply this with M — L\. Thus 

£ = L 

In this sum, we group the L G £ such that L + L is equal to a given 
Mi G £ together. 

First we consider the case when dimLi — dimL = 1. Then Mi is either 
Lq or Li, hence we obtain 

/(L ,L ) + /(L ,L 1 ) = 1. 

Since f(L ,L ) = 1 by (a), we obtain f(L ,Li) = as required. 
Next we consider the case when dimLi — dimLo > 1. We obtain 

£ /(Lo,Mi) = l. 

{M16A LoCMiCLi} 

For Mi = L , we have f(L ,L ) = 1 by fa/ For Mi S Li, we have 
/(L ,Mi) = by induction hypothesis. Hence there remains only the term 
f(L , Li) which must be equal to 0. □ 

We study now the Euler-Maclaurin expansion of a Barvinok valuation. 

Theorem 17 Let p C V be a rational polyhedron and let f be a face of p. 
Let £ be a finite family of rational vector subspaces of V , stable under sum. 
Let p be a patchwork function on £, and let S^' p ' = ^2 LeC p(L)S L . 
Assume that lin(f) belongs to £. Then 

^)(t(p,f))=/iW(t(p,f)). 

In other words, the j-term in the local Euler-Maclaurin expansion of S^' p \p) 
coincides with that of S(p). 

Proof. In the sum ^Le£P(-^)A iL+lm( ' f ^ lm< '^(t(P) f))> we group the terms for 
which L + lin(f) is equal to a given L\ together. By Lemma [16] we obtain 
/i {0} (t(p, f)) for Li = lin(f) and otherwise. □ 
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Corollary 18 Let p C V be a rational polyhedron. Let < k < d. Let C 

be a finite family of rational vector subspaces of V , stable under sum, such 
that lin(f) 6 C for every k-dimensional face f of p and let p be a patchwork 
function on C. 

• Let < k < d. Then the meromorphic function 

s(p)(0-s^\p)(0 

has lowest degree > —k + 1. 

• Let k = 0. Then 

s(p)(0 = s^\p)(0- 

Proof. By Theorem [TTJ, the local Euler-Maclaurin expansion of the dif- 
ference involves only faces of dimension < k. 

s(p)(0-s^\p)(0 = Yl (/i {0} (t( P ,f))(O -^' p) (t(P,f))(0)/(f)(0. 

{fe^(p), dimf<fc} 

For a face of dimension j, the function /(f) (£) is homogeneous of degree —j. 
Multiplied by the holomorphic function n^(t(p, f))(£) — p c {t{p, f))(£)> i ts 
lowest degree can only increase. □ 

Remark that the statement of Corollary [TS] above does not involve any 
scalar product. In the next section, we will show that this corollary implies 
our main Theorem (Theorem |20{) . 

4 Ehrhart quasi-polynomial 

Let p be a rational polytope and let h(x) be a polynomial function on V. Let 

s(p,h)= h w- 

x&pnz d 

When p is dilated by a non negative integer n, we obtain the Ehrhart quasi- 
polynomial of the pair (p, h) 

d+N 

S(np,h) = J2E m (p,h,n)n m , 

m=0 
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where N = degh. The coefficients E m (p,h,n) are periodic functions of 
nGN, with period the smallest integer q such that qp is a lattice polytope. 

If an integer r < d is fixed, and h = 1, Barvinok [2 J proved that the r + 1 
highest Ehrhart coefficients Ed(p, 1, n), . . . , _Ed_ r (p, 1, n) of 5 (rap, 1) can be 
computed in polynomial time with respect to d, when p is a rational simplex. 

Let L C V be a rational vector subspace. Denote the projected lattice 
on V/L by A^. Consider the mixed valuation 

S L (p,h) = J2 [ h(x)dx. 

As shown by Barvinok, and as we will see here, we can use linear combination 
of these mixed valuations to approximate S(np, h) when n is big. 

For any polyhedron a, we define the meromorphic function S L (a, h)(£) G 
A4i(V*) similarly to S L (a,h). For ( G such that the sum converges, we 
have 

S L (a,h)(g) = [ h(x)e^ x) dm L (x). 

Remark 19 It is clear that S L (a, h)(£) = h(d^) ■ S L (a)(0- 

If p is a polytope, then S L (p, h)(£) is regular at £ = and S L (p, h)(0) = 
S L (P,h). 

For a family C and a function L \— > p(L), L G C, we define 

&^\a,h)® = 'E t p{L)S L (a,h)it). 

Lec 

If p is a polytope, and we dilate by n G N, we have again a quasi- 
polynomial 

d+N 

S^\np,h)(0) = J2E m (C,p,p, h,n)n m . 

We can replace the quasi-polynomial S^ c ' p \np, h)(0) by q legal polynomi- 
als in the variable u, by splitting N into classes modulo q. Writing n = qu + k, 
for k = 0, . . . , q — 1, we obtain the polynomial function of u: 

d+N 

S^\{qu + k)p, h)(0) = £ *# } (4 Pi Pi h > m - 

m=0 
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We briefly recall how the usual Ehrhart quasi-polynomial of a polytope 
can be computed using Brion's theorem. We will then use a similar method in 
order to compute efficiently the r + 1 highest coefficients only, using Barvinok 
valuations. 

Let V(p) be the set of vertices of p. For each vertex s, let c s be the cone 
of feasible directions of p at s, so that the supporting cone at s is s + c s . By 
Brion's theorem [5], we have 

S(p,h)(t)= E S(a + c,h)(£)- 

s G V(p) 

Let n G N and consider the dilated polytope np. The supporting cone at the 
vertex ns is ns + c s . Let q G N such that qp is a lattice polytope and fix 
k &N, < k < q — 1. Let n = qu + k. As gus is an integral point, we have 

S((qu + k)s + c„ h)(0 = e qu ^ s) S c ^(ks + c„ h)(£). 

Expanding in powers of u, we obtain 

S((qu + k)p,h)(0 = Ys Umq -n E (Z,s) m S(k8 + c„h)(Z). 

z — ' ml z — ' 

It follows that for each m, the sum of meromorphic functions 

q — (Z,s) m S(ks + c s ,h)(0 
ml L — ' 

sev(p) 

is actually analytic. Its value at £ = is obtained by taking the zero degree 
term. We obtain 

S((qu + k)p,h)(0) = J2 E ™^V « m > 

m>0 

with 

£S ) (P^) = ^ E (^s) m S(ks + c s ,h) hm] (0. 

The right-hand side of this relation, a priori a meromorphic function of £, is 
actually constant. Moreover, we have S^/cs + c s , h)[- m ](£) = if m>d + N, 
hence the Ehrhart quasi-polynomial has degree < d + N. 
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We apply now Brion's theorem to S c,p (p, h)(£). We obtain 

S £ "(p,h)(Z)= E s c,p (s + c s ,h)(0. 

sev(p) 

For reasons to be explained later on, instead of one family C, we take a family 
of subspaces C s for each vertex s. Let p s : C s — > Z be a function on £ s . We 
denote now by (£,p) the map s h- > (C s ,p s ). 
We define: 

(9) ^(p,/ i )(0= ]T fi^*( a + c.,/0(O. 

If the family does not depend on s, (C s , p s ) = (£ 0) Po) f° r every vertex s, 
then, by Brion's theorem, we have 

^(p,/ i )(0 = ^°(p,/ i )(0- 

We dilate ©. Let n = qu + k. We obtain 

B^((gu + A;)p ) /i)(0= E e^S^'CAw + c s , h){t). 

sev(p) 

Expanding in powers of u, we obtain 

(10) B L >o((qu + k)p, h)(0 = J2 u m E%\C, p, p, h)(0 

m>0 

with 

(11) E%\£.,p,p,h){0 = —< E ^s) m S c ^(ks + c s ,h)(0- 

mi z — ' 

If the family does not depend on s, (C a , p s ) = (Co, po) for all vertices, then 
B c ' p (p, h)(£) = S c °' P0 (p,h)(£) is analytic near £ = 0, and so are the coeffi- 
cients ffTTj) . 

On the contrary, if we take a different family C s for each vertex s, the 
coefficient Em(C, p, p, /i)(£) of u m in ([TO]) is no longer analytic near £ = 0, in 
general. However, the meromorphic function £ i— > Em\jC, p, p, /&)(£) belongs 
to .M^V*), thus it has a term of degree with respect to £, given by 

(12) ^)(£,p,p,/ i ) [0] (0 = ^y E (ts) m S C ^(ks + C s ,h) hm] (0. 

' sev(p) 
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For a family (C,p) as described in the next theorem, it turns out that, for 
large m, this zero-degree part Em (A p, p, /i)[o](£) is actually analytic, hence 
constant, and its value is equal to the m-th Ehrhart coefficient Em\p, h) of 
S((k + qu)p,h)(0). 

Theorem 20 Let p be a rational polytope in a rational vector space of di- 
mension d. For each vertex s of the polytope p, let c s be the cone of feasible 
directions ofp at s, so that the supporting cone at s is s + c s . For each vertex 
s, let C s be a finite family of rational vector subspaces of V , stable under 
sum, such that lin(f) belongs to C s for every face f of codimension r of the 
cone c s , and let p s be a patchwork function on C s . Let q G N such that qp is 
a lattice polytope and fixk&N,0<k<q — 1. Let h(x) be a homogeneous 
polynomial of total degree N. 

Then, form > d + N — r, the zero-degree term Em\c,p, p, fo)[o](0 defined 
by (TTjj) is regular near £ = 0, hence constant. Its value is the coefficient 
Em\p,h) of u m in the Ehrhart quasi-polynomial 

d+N 

S((k + qu)p, h)(0) = h ^ = E ^-(P, h). 

xe((k+qu)p) nA rn=0 

Proof. We first consider the case h(x) = 1. We have, for every m > 0, 

a m 

E{ ™ ) (P' 1 ) = h E (^s) m S(ks + c s ) hm] (0 

lit. 

sev(p) 

where the right-hand side is actually a constant function of £. For m > 
d — r — 1, we have, by Corollary [T8| 

S(k S + C s )(0[-m] = S Ca ' P °(ks + C s )[_ m ](0- 

This proves the theorem when h(x) = 1. The case of a non constant 
polynomial h(x) is quite similar. If h(x) = x^ 1 . . . x^ d , we just have to 
replace the meromorphic functions S(ks + c s )(£) and S Cs ' ps (ks + c s )(£) by 
their derivatives under d^ 1 . . . d f d . □ 

si 

If for each vertex s, we take C s = £, the full collection generated by all 
r codimensional faces of p, we obtain Corollary [2T1 below, that is Barvinok's 
theorem [2J, with an extension to the sum of values of any polynomial h(x) 
over the set of integral points of a rational polytope (Barvinok considers only 
the case h(x) = 1). 
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Corollary 21 Let p C V be a rational polytope and let h(x) be a polynomial 
function on V. Let £ be a finite family of rational vector subspaces of V , 
stable under sum. Assume that lin(f ) belongs to £ for every face f of codimen- 
sion rofp. Let p be a patchwork function on £ and let S c ' p = ^2 LeC p(L)S L . 
Then the r + 1 highest Ehrhart coefficients of S(tp,h)(0) and S c ' p (tp, h)(0) 
are equal. 

The point in taking a family £ s which depends on the vertex s lies in the 
case where p is simplicial. In this case, we can take £ s to be just the set 
of subspaces lin(f), for all faces f of codimension < r of the supporting cone 
c s at vertex s. This family is stable under sum. Moreover the patchwork 
function on £ s is simple, (Lemma [T5j) and the computation of the function 
S L (ks + c s )(£), when L e £ s , is immediately reduced (Example [6]) to the 
computation of a function S(a)(C,) for a simplicial cone a in a rational vector 
space of dimension smaller or equal than r. When p is a simplex, we obtain 
in this way a method for computing the r + 1 highest Ehrhart coefficients for 
the pair (p, h). 

5 Local Euler-Maclaurin formula for mixed 
sums 

Finally in this last section, we discuss an application of the existence of the 
coefficients p L (Theorem [8]) in the line of [1]. 

Let p be a rational polytope in a rational vector space V of dimension d 
and let h{x) be a polynomial function on V. Let L be a rational subspace of 
V. Consider the mixed sum 



As in [1], we associate to the analytic function /i L (t(p,f)) a constant 
coefficients differential operator (of infinite order) on V. 

Definition 22 Let f be a face of p. We denote by D L (p,f) the differential 
operator on V associated to analytic function fi L (t(p, f)): 




D L (p,m) ■ = p L (t(p,f)))(Oe^ x) . 
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The operators D L (p, f) are local, that is they depend only of the trans- 
verse cone t(p, f) of p along f, and they involve only derivatives in directions 
orthogonal to the face f. We can state the following theorem with the same 
proof as in [3]. 

Theorem 23 (Local Euler-Maclaurin formula) Let p be a polytope in V . For 
any polynomial function h(x) on V , we have 



where the integral on the face f is taken with respect to the Lebesgue measure 
on < f > defined by the lattice A D lin(f) . 
In particular, for h = 1, we obtain 



Let us dilate the polytope p by a non negative integer n. If f is a face of 
p, let qf be the smallest positive integer such that ?f < f > contains integral 
points. Define D(p,f, n) = D(np,nf), if n > 0, and D(p,f, 0) = D(qft,qtf). 
The function n \— > D(p, f, n) is periodic of period qf. 

Proposition 24 Let p be a rational polytope and h a polynomial function of 
degree N on V. Then, for any integer n > 0, we have 



(13) 




(14) 




(15) 




Furthermore, if f G .F(p), have 




dim f+JV 



where the coefficients Ei(p, h, f, n) are periodic with period qf. 
Hence the Ehrhart coefficients are given by 



Ei(p,h,n)= E m (p,h,f,n). 
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When we apply the last proposition to the function h(x) = 1, we obtain 
(16) S L (np,l)= /i L «P,f))(0)vol(f)n dim f. 

As /z L (nt(p, f)) is invariant by integral translations, the function fj, L {nt(p, f))(0) 
is of period q$. 
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